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ABSTRACT 

A standard binary microlensing event lightcurve allows just two parameters of the lensing 
system to be measured: the mass ratio of the companion to its host, and the projected sep- 
aration of the components in units of the Einstein radius. However, other exotic effects can 
provide more information about the lensing system. Orbital motion in the lens is one such 
effect, which if detected, can be used to constrain the physical properties of the lens. To 
determine the fraction of binary lens lightcurves affected by orbital motion (the detection 
efficiency) we simulate lightcurves of orbiting binary star and star-planet (planetary) lenses 
and simulate the continuous, high-cadence photometric monitoring that will be conducted by 
the next generation of microlensing surveys that are beginning to enter operation. The effect 
of orbital motion is measured by fitting simulated lightcurve data with standard static binary 
microlensing models; lightcurves that are poorly fit by these models are considered to be de- 
tections of orbital motion. We correct for systematic false positive detections by also fitting 
the lightcurves of static binary lenses. For a continuous monitoring survey without intensive 
follow-up of high magnification events, we find the orbital motion detection efficiency for 
planetary events with caustic crossings to be 0.061 ± 0.010, consistent with observational re- 
sults, and 0.0130 ± 0.0055 for events without caustic crossings (smooth events). Similarly for 
stellar binaries, the orbital motion detection efficiency is 0.098 ±0.011 for events with caustic 
crossings and is 0.048 ± 0.006 for smooth events. These result in combined (caustic crossing 
and smooth) orbital motion detection efficiencies of 0.029 ± 0.005 for planetary lenses and 
0.070 ±0.006 for stellar binary lenses. We also investigate how various microlensing parame- 
ters affect the orbital motion detectability. We find that the orbital motion detection efficiency 
increases as the binary mass ratio and event time-scale increase, and as impact parameter 
and lens distance decrease. For planetary caustic crossing events, the detection efficiency is 
highest at relatively large values of semimajor axis ~ 4 AU, due to the large size of the res- 
onant caustic at this orbital separation. Effects due to the orbital inclination are small and 
appear to only significantly affect smooth stellar binary events. We find that, as suggested by 
Gaudi (2009), many of the events that show orbital motion can be classified into one of two 
classes. The first class, separational events, typically show large effects due to subtle changes 
in resonant caustics, caused by changes in the projected binary separation. The second class, 
rotational events, typically show much smaller effects, which are due to the magnification 
patterns of close lenses exhibiting large changes in angular orientation over the course of an 
event; these changes typically cause only subtle changes to the lightcurve. 

Key words: gravitational lensing - binaries: general - stars: low mass, brown dwarfs - plan- 
etary systems - Galaxy: general 



1 INTRODUCTION 

The current gravitational microlensing surveys, OGLE (Udalski 
|2003| > and MOA < |Hearnshaw et al.|2005) discover ~ 700 unique 
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microlensing events per year, of which, of order ten percent show 
signatures of lens binarity. A small fraction of these, those with a 
high probability of planet detection, are followed-up by a number 
of follow-up teams, which intensively monitor the events for the 
signatures of planets. In the coming years this strategy will be aug- 
mented and extended by a strategy of continuous, high cadence sur- 
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veys, performed by a global network of wide field telescopes. Such 
a network will monitor all the microlensing events it discovers with 
a cadence similar to that achieved by the follow-up networks for a 
handful of events today. 

The lightcurve of a standard static binary lens, in which the 
lens components are fixed and the source follows a straight path, 
can be described by a minimum of seven parameters. Only three of 
these parameters contain physical information about the lens sys- 
tem. Two are dimensionless parameters: the mass ratio q, and the 
projected separation of the lens components d, measured in units of 
the Einstein radius. The third, the Einstein time-scale of the event 
tE, is the time taken for the source to cross one Einstein radius 



t E = 



(1) 



where v t is the relative projected lens-source velocity and te is the 
Einstein radius. This is defined as 



te 



4G 



x(l - x)D s M, 



(2) 



where x = Di/D s is the ratio of the lens distance D\ to the source 
distance D s and M is the total mass of the binary. Of the other four 
parameters, three are purely geometrical, and the final parameter is 
the unlensed source flux. 

The mass ratio and separation are closely related to the most 
interesting properties of the binary, the component masses and the 
semimajor axis of the orbit. They can be measured very accu- 
rately from a lightcurve, but only describe the binary's properties 
in terms of ratios relative to the typical physical scales of the sys- 
tem. The Einstein crossing time-scale £e contains information on 
these scales, but this information is wrapped up in a three-fold de- 
generacy (the so-called microlensing degeneracy) between the total 
binary mass, the lens distance and the source velocity. It is also de- 
pendent on the source distance, but this is usually well constrained 
by measurements of the baseline flux. To gain any more knowledge 
of the lens system requires that this degeneracy be broken, either by 
the detection of higher order effects in the event lightcurve, or by 
detection of the lens flux and proper motion as the lens and source 
separate]^] These detections yield measurements of the lens dis- 
tance and source velocity respectively, allowing the lens mass to be 
solved for JGould & Loeb|1992| [Bennett et al.|2006l>. H igher order 
effects, such as finite source effects ( |Witt & Mao|1994l|Nemiroff &| 
|Wickr amasinghe 1994||Alcock et al.|1997|> and m icrolensing par- 
allax ( |Refsdal||1966||Gould||1992||Alcock et al.|[l99"5} , allow the 
microlensing degeneracy to be broken or reduced through measure- 
ment or constraint of some of the parameters that are combined in 
tE - For example detections of finite source effects and microlensing 
parallax in the same event yield two independent measurements of 
the angular Einstein radius 6e = te/D\, which allow the source 
velocity and lens distance to be eliminated, and the lens mass de- 
termined (e.g. |An & Gould|2001| >. 

Orbital motion of the binary lens is another such higher order 
effect. If the binary lens components are gravitationally bound, they 
will orbit each other, and their projected orientation will change as 
a microlensing event progresses. As the magnification pattern pro- 
duced by a binary lens is not rotationally symmetric, the change in 
orientation may be detectable in the lightcurve of the event. If the 
orbit is inclined relative to the line of sight, then the projected sep- 
aration of the lens components will also evolve, causing changes in 



1 Throughout we will use the terms lens motion and source motion inter- 
changeably. 



the structure of the magnification pattern, which again may be de- 
tectable. In a small fraction of binary microlensing events we can 
expect to see the effects of this orbital motion in their lightcurves, 
though this is the first work that attempts to quantify this fraction. 
If orbital motion can be detected in a microlens it can provide con- 
straints on the mass of the lens, and information about the binary 
orbit. 

To date, six binary microlensing events have shown strong ev- 
idence of orbital motion in the lens system. The first, MACHO-97- 
BLG-41 was a stellar mass binary. Modelling of the event was only 
able to measure the change in the projected angle and separation 
of the binary in the time between two caustic encounters, but was 
unable to constrain the orbital parameters ( Albrow et al. 2000}. The 
second event, EROS-BLG-2000-5, had very good lightcurve cov- 
erage, which allowed the measurement of the rates of change of 
the binary's projected separation and angle; these measurements 
were then used to obtain a lower limit of the orbit's semimajor axis 
and an upper limit on the combined effect of inclination and ec- 
centricity ( |An et al.|2002| >. The third and fourth examples, OGLE- 
2003-BLG-267 and OGLE-2003-BLG-291 both seem to show or- 
bital motion effects ( Jaroszynski et al. 2005 ). However, only OGLE 
survey data was used in their analysis, without follow-up measure- 
ments, so the lightcurve coverage was not ideal. Combined with 
parallax measurement, the masses of both binary lenses were con- 
strained, but no constraints could be placed on the orbits ( Ja roszyn^| 
|ski et al.|2005| l. In each of these four cases, the ratio of the com- 
ponent masses is large (near unity), indicative of the lens systems 
being binary stars, however, orbital motion has recently been mea- 
sured in two events involving planetary mass secondaries. After the 
paper was submitted, two further events have been shown to dis- 



play orbital motion effects: OGLE-2005-BLG-153 (Hwang et al. 
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and OGLE-2009-BLG-092/MOA-2009-BLG-137 ( |Ryuetal. 



OGLE-2006-BLG-109 was an event involving a triple lens, 
with analogues of Jupiter and Saturn orbiting a ~ 0.5Mq 
star ( |Gaudi et al.[2 008 1. The lightcurve of the event had extremely 
good coverage, and showed multiple features, allowing the orbital 
motion of the Saturn analogue to be detected. The detection was 
so strong that the semimajor axes of both planets could be strongly 
constrained ( Gaudi et al. 2008 1. A more complete analysis of the 
event, incorporating measurements of the lens flux and orbital sta- 
bility constraints, carried out by Bennett et al. (2010), tightly con- 
strained four out of six Keplerian orbital parameters of the Saturn 
analogue, and weakly constrained a fifth. The planet OGLE-2005- 
BLG-071Lb is a ~ 4 Jupiter mass planet orbiting a ~ 0.5Mq 
star ( |Udalski et a l. 2005). Measurements of the orbital motion 
in this event have allowed some constraints to be placed on the 
planet's orbit (Dong et al. 2009). In all six events other higher or- 
der effects have also been detected, most notably microlens paral- 
lax and finite source effects, which are detected in all the events, 
and in each case allow the measurement of the lens mass. 

Despite these detections, there has been relatively little the- 
oretical work on orbital motion in microlensing, likely due to the 
traditional assumption that the effects of orbital motion on a binary 
microlens lightcurve will be small and in most cases negligible (e.g. 
|Mao & Paczyris ki|1991| >. The problem was first considered in de- 
tail by Dominik ( 1998 1, who concluded that in most microlensing 
events the effects of lens orbital motion were likely to be small, 
though in some cases lightcurves could be dramatically different. 
Dominik ( 1998) points out that the effect is most likely to be seen 
in long duration binary microlensing events with small projected 
binary separations. |Ioka, Nishi & Kan-Ya| |l999) also studied the 
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problem, and pointed out that the effect of binary lens rotation 
is likely to be important in self-lensing events in the Magellanic 
clouds. Rattenbury et al. (20021 showed that orbital motion could 
affect the planetary signatures seen in high-magnification events. 

The six microlensing events that display orbital motion make 
up a significant fraction of the few tens of binary microlensing 
events that have been modelled (e.g. jAlcock et al.|2000||Jaroszyn-| 
[ski]|2002t [Jaroszynski et al.||2004| |2006| |Skowron et aLpOOTf T 



which begins to shed doubt on the previous conclusion that lens or- 
bital motion is likely to be unimportant in most binary events. The 
two planetary events constitute approximately 15 percent of the en- 
tire published microlensing planet population. These observations 
motivate us to revisit the question: how likely are we to see lens or- 
bital motion in a microlensing event? This question is made espe- 
cially pertinent in the context of the next generation of high cadence 
microlensing surveys which will make the exquisite lightcurve cov- 
erage of EROS-BLG-2000-5 and OGLE-2006-BLG-109 the norm 
rather than the exception. To gain a better understanding of how fre- 
quently orbital motion affects microlensing lightcurves we simulate 
a large number of microlensing events caused by orbiting binary 
lenses. We also investigate the factors that affect this frequency. 

The structure of the work is as follows. In Section [2] we will 
review the basic theory of binary microlensing, and the effects of 
orbital motion on such lensing systems. Section [3] describes our 
simulations of microlensing events and Section Bfdescribes how 
we measure the effects of orbital motion. In Section [5] we present 
the results of the simulations. We discuss the results in Section [6] 
and conclude in Section|7] 



2 MICROLENSING WITH ORBITING BINARIES 
2.1 Binary microlensing 

The lens equation of a binary point mass gravitational lens de- 
scribes the mapping of light rays from the source plane to the image 
plane, and can be written in complex form (Witt 19901 

1/1 



1 + q \ z — z\ 



+ 



Z — Z2 



(3) 



where z s = x s + iy s is the complex coordinate in the source plane, 
z — x + iy is the complex coordination the image plane, q ^ 1 is 
the mass ratio of the secondary mass to the primary, z\ = x\ + iyi 
and Z2 = X2 + iy2 are the complex coordinates of the primary 
and secondary lens respectively, and bars represent complex con- 
jugation. All lengths have been normalized to the Einstein radius 
of the total lensing mass. The positions of the images z for a given 
source position are found by solving this equation, which can be 
rearranged into a fifth-order complex polynomial in z. The total 
magnification of the images is given by the ratio of image areas to 
the source area. This information is contained in the Jacobian of the 
lens mapping J, and the magnification is given by 

^detJ' (4) 
where det J is the determinant of the Jacobian and is given by 

2 

a \ 

+ 



det J =1 



1 + q V(^-2i) 2 (z-zi) 2 



(5) 



2 The symbol x used here should not be confused with that representing 
the ratio of lens to source distances. Its meaning should be clear from the 
context in which it is used, and it will not be used again in this context 
without a subscript. 
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Figure 1. Plot showing the d-q plane separated into three regions where the 
caustics take on close, resonant and wide topologies, with increasing d; the 
dashed lines, d c (q) and d w (q), separate the regions of different topology. 
The solid lines are examples of each of the caustic topologies, drawn to the 
same scale, for a binary lens with q = 0.1, and d = 0.7 (close), d = 1.05 
(resonant) and d = 1.75 (wide). The filled circles show the positions of 
the lenses for each topology, with the primary (more massive) lens being 
positioned leftmost in all cases. After Cassan 1 2008 1. 



It is possible for det J to be zero; if this occurs, the magnifi- 
cation of a point source will be infinite at the point in the source 
plane where det J = 0. This occurs when the quantity within the 
modulus sign in equation (|5j lies on the unit circle in the Argand 
diagram. Any point in the image plane that obeys det J = is a 
critical point. The set of critical points form a set of closed curves 
called critical curves. The critical curves can be found by solving 
the critical curve equation 



(z — Zi) (z — z 2 y 



= (1 + q)e 



(6) 



where is a parameter, such that when swept over ^ (f> < 2n, the 
solutions z draw out the critical curves. These can then be mapped 
through the lens equation {3} onto the source plane to form caustics. 
A point source that lies on such a caustic will be infinitely mag- 
nified, but this unphysical magnification remains finite for physi- 
cal, finite sources. The caustics are characteristically made up of 
smooth curves, fold caustics, that meet at cusps. 

The number and shape of caustics is determined by just two 
lens parameters, the mass ratio q and the projected lens separa- 
tion d in units of Einstein radii. There are three possible topologies 
that the caustics can assume: close, resonant and wide; examples 
of each are shown in Figure [T] The lines that delimit the different 
topologies in the (d, q) plane, also shown in the figure, are given 
by |Erdl & Schneider| (1991) 1: 



il' 



(1 + g) 
27q 



<l-dt) 3 



which separates regions of close and resonant topology, and 

w ~ 1+g ' 
which separates regions of resonant and wide topology. 



(7) 



(8) 
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2.2 Orbital motion in a binary microlens 



3.1 The Galactic model 



The lightcurve of a microlensing event can be considered as a one- 
dimensional probe, by the source, of the two-dimensional magnifi- 
cation pattern produced by the lens. The magnification pattern of 
a single lens is rotationally symmetric about the position of the 
lens, but the magnification pattern of a binary lens is more compli- 
cated, containing strong caustic structures that exhibit a refiectional 
symmetry about the binary axis, the axis connecting the lens com- 
ponents. However, far away from the caustics, the magnification 
pattern can resemble that of a single lens. 

As the lens components orbit each other, their position an- 
gle and their projected separation can change. These changes cause 
changes in the orientation and structure of the magnification pattern 
respectively. It is clear, however, that only if the source traverses re- 
gions of the magnification pattern that differ significantly from that 
of a single lens, will it be possible to detect these effects of orbital 
motion. For the effects to be measurable the lightcurve of the event 
must be affected in a significant way, that is not reproducible by 
a static binary lens model. It is also possible to detect the effect 
of orbital motion by showing that a static model is less physically 
plausible than an orbiting model, but this will usually require fur- 
ther information about the event, such as an independent constraint 
on the lens mass. 

The effects of orbital motion on a lightcurve can also be mim- 
icked by other higher order effects, especially parallax and xallarap. 
Parallax effects are caused by the motion of the earth about the sun, 
and cause the source to take an apparently curved path through the 
magnification pattern (e.g. Smith, Mao & Paczyhski 20031. In the 
case of xallarap, the source travels along a curved path through 
the magnification pattern as a result of binary orbital motion in the 
source system (|Griest & Hu|1992|[Paczynski|19971|Dominik[T998l 
|Rahvar & D ominik 2009). These curved paths can look very simi- 
lar to those taken by the source in the rotating binary lens centre of 
mass frame, and hence it can sometimes be difficult to identify the 
true cause of the effect. 



3 SIMULATING A HIGH CADENCE MICROLENSING 
SURVEY 

The major aims of this study are two-fold: firstly to determine the 
fraction of microlensing events that will be affected by orbital mo- 
tion, as seen by the next generation microlensing surveys; and sec- 
ondly, to investigate the factors that affect the detectability of or- 
bital motion, to aid the targeting of such events without resorting to 
exhaustive modelling efforts. To achieve the first goal, the various 
factors that go into the observation of a microlensing event should 
be simulated, accurately modelling the observing setup, the distri- 
butions of planetary and binary star lens systems, and the distribu- 
tion of the sources and lenses throughout the Galaxy. To achieve 
the second goal we must simplify the parameter space we investi- 
gate as far as possible, without removing essential elements from 
the model, so as to allow a clear interpretation of the results. 

To balance these somewhat contradictory requirements we 
choose to accurately simulate ideal photometry and use a semi- 
realistic model of the Galaxy, while investigating a logarithmic dis- 
tribution of companion masses and separations. This allows us to 
use our simulations to gain a good order of magnitude estimate of 
the results expected from future surveys, whilst simultaneously in- 
vestigating the factors that have the largest impact on the detection 
of orbital motion over a relatively uniform parameter space. 



To simulate the kinematic and distance distributions of the source 
and lens populations we assume a simplistic bulge and disk model 
of the Galaxy. We assume the source to be located in the bulge, at 
a fixed distance D s = Rg — 8 kpc, in the direction of Baade's 
Window, where Ro is the distance to the Galactic centre. The lens 
distances are distributed according to the stellar density distribution 
of Model II of Binney & Tremaine ( 2008 1, which consists of a thin 
and a thick exponential disk and an oblate spheroidal bulge with 
a truncated power-law density distribution. The kinematics of our 
Galactic model are based on that of |Han & Gou ld ( 1995bj), who 
describe the kinematics of a stellar disk and a barred bulge. The 
distribution of relative source velocities dn/d« t is dependent on 
the transverse velocities of the lens, source and observer, and their 
corresponding velocity distributions. The observer is assumed to 
follow the Galactic rotation at the position of the Sun, and there- 
fore has a velocity (i>o/, «o,i>) = (225.2,7.2) km s~ x in the di- 
rections of Galactic coordinates (£, b), once the Solar apex motion 
is included. The source and lens are assumed to follow the Galac- 
tic rotation, with an additional random component. Their velocities 
have the form, in the directions £ and b, 



Vl — l>rot + l>rand,f, Vb — t>rand,b, 



(9) 



where « ro t is the rotational component of the velocity, and v ran a,i 
and Urand.b are random velocities in the directions I and b respec- 
tively. The rotation curve of the bulge is assumed to be flat beyond a 
distance of 1 kpc from the Galactic centre, and that of a solid body 
within 1 kpc. Therefore, the rotational velocity component v TO t for 
bulge stars is 



Wrot 



ifi?<lkpc 
x if R ^ 1 kpc 



(10) 



where v ma x = 100 km s 1 is the maximum rotational velocity 
of the bulge, and R = V ' X 2 + Y 2 and (X, Y, Z) is a galac- 
tocentric coordinate system with the X-axis increasing towards 
the observer and the Z-axis pointing out of the Galactic plane; 
for the disk v vot = 200 km s~ . The random velocity compo- 
nents are assumed to follow Gaussian distributions, with disper- 
sions taken from |Han & Gould| (1995a| >. These dispersions are 
(o"£, ctj,) = (30, 20) km s~ for the disk and (ax, ay, oz) = 
(110, 82.5, 66.3) km s _1 for the bulge. From these quantities, the 
relative transverse velocity of the source v t (the quantity we are in- 
terested in) can be calculated from the relative velocities in the £ 
and b directions vi and Vb as 



v t = 



where (e.g. |Han & Gould|1995b) 



2 + <, 



vi,b = (v\ — v )t,b + x(v — v a )t,b, 



(11) 



(12) 



and v , vi and v s are the observer, lens and source velocities re- 
spectively, in the directions £ and b. 

The final distribution of lens distances and velocities takes into 
account the dependence of the event rate T oc -\/x(l — x)v t on the 
distribution of each parameter. While the kinematic and density dis- 
tributions are produced from different Galactic models, they quali- 
tatively reproduce the observed Einstein crossing time-scale distri- 
bution, shown in Figure|2] including its asymptotic behaviour ( Mao 
|& Paczynski|1996| >. 
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Figure 2. The Einstein time-scale distribution for ~ 50000 simulated 
events. The solid line and data points show the simulated data, and the 
dashed lines show lines of slope 3 and —3, the expected asymptotic be- 
haviour of the distribution. 



3.2 The microlensing events 

When observing a microlensing event, it is often the case that the 
light of the source being magnified is blended with that of nearby 
stars in the field. The amount of blending can be quantified by a 
blending fraction / s , which we define to be the fraction of the total 
flux of the observed blend that the source contributes when unmag- 
nified, such that the time dependent magnitude of the blend is 



7(f) =/ b -2.51og(/ s A(t) + (1 -/.)), 



(13) 



where 1^ is the baseline magnitude of the observed blend when the 
source is unmagnified, and A(t) is the magnification caused by the 
lens. 

The distribution of baseline magnitudes and blending frac- 
tions is drawn from simulations of blending effects by Smi th et ah] 
( 2007 1, who perform photometry on mock images of typical Galac- 
tic bulge fields with high stellar density. Specifically we calcu- 
late the blending fraction and baseline magnitude for each event 
from the input and output magnitudes of source stars drawn from 
their simulation with 1.05 arcsec seeing and input stellar density 
of 133.1 stars arcmin -2 , before any detection efficiency cuts are 
made to the catalogue. As the phenomenon of negative blending, 
the source apparently contributing a fraction f s > 1 to the total flux 
of the blend ( |Park et al.|2004||Smifh et al.|2007) , is poorly under- 
stood, we only include sources with moderate negative blending, 
requiring that f s < 1.2. 

The mock images are produced by Smith et al. ( 2007 1 using 
the method of |Sumi et al.| |2006 l, drawing stars from the Hubble 
Space Telescope /-band luminosity function of Holtzman et al. 
( 1998 1, adjusted to account for denser fields and brighter stars us- 
ing OGLE data. Extinction was accounted for using the extinction 
maps of Sumi (20041, and the baseline magnitudes were measured 
using the standard OGLE pipeline based on DOPHOT (Schcchtcr 
1993| >. Further details can be found in section 3 of the |Smith 



et al. 



et al. 



(2007 ) paper, and references therein. 



The lens systems are composed of a primary of mass M± , and 
secondary of mass A/2. The primary's mass is drawn from a broken 
power-law distribution 

dn M (a+0.B) _ f -L3 Ml SC Mbreak 

dMx 1 '""I -2.0 Mi > M brcak 
with lower and upper limits of 0.05Mq and 1.2Mq respectively, 



(14) 



and where M break = 0.5Mq. The addition of 0.5 to the power-law 
index is to account for the dependence of the microlensing event 
rate on the mass of the lens. We do not include a population of stel- 
lar remnant lenses, such as white dwarfs, neutron stars and black 
holes. The mass ratio q of the secondary to the primary is drawn 
from a logarithmic distribution, with limits 10~ 2 ^ q < 1 for stel- 
lar binary lenses and 10~ 5 ^ q < 10~ 2 for planetary lenses. Note 
that for lower mass primaries, the distribution of stellar binary mass 
ratios does include secondaries with masses as low as ~ 5Mj up iter, 
well into the planetary mass regime, and the lower limit of the plan- 
etary mass ratio distribution implies a secondary of ~ for a 
0.3Mq primary. 

The components of the lens orbit their combined centre of 
mass in Keplerian orbits, of semimajor axis a, distributed loga- 
rithmically over the range 0.1 AU ^ a < 20AU. These orbits are 
inclined to the line of sight, with inclination angles distributed uni- 
formly over a sphere. For binary stars we performed two sets of 
simulations, one with zero eccentricity e, and another with bound, 
eccentric orbits with eccentricities distributed uniformly over ^ 
e < 1. 

The source trajectories were parametrized by the angle of the 
source trajectory relative to the binary axis 6a, at the time of closest 
approach to, and the impact parameter uq, the projected source-lens 
separation in units of Einstein radii at to- We set to = 0, for sim- 
plicity, and 9q and no were distributed uniformly over the ranges 
^ do < 2-7T and —1.5 ^ ito < 1.5 respectively. 



3.3 Simulation of photometry 

In the hunt for planets, the proposed next generation of microlens- 
ing survey will consist of a (potentially homogeneous) network of 
telescopes located throughout the southern hemisphere such that 
the target fields in the Galactic bulge can be monitored contin- 
uously during the times when the bulge is observable. The tele- 
scopes will have diameters between 1.3-2.0 m and fields of view 
1.4-4.0 square degrees. They will operate at a cadence of approx- 
imately 10 minutes, and are expected to discover several thousand 
microlensing events per year. An example is KMTNet, a network of 
three identical 1.6m telescopes due to enter operation in 20 1 4 ( Kim 
|et al.|20i0| >. Such surveys can operate effectively without the need 
for intensive follow-up observations due to their high cadence and 
continuous coverage. However, it is likely that the survey/follow-up 
observing paradigm will persist, with low cadence surveys monitor- 
ing far larger areas of sky. 

Unfortunately the effects of the weather, amongst other 
things, makes completely continuous, high-cadence observations 
unachievable in reality. Rather than including complicated models 
of these effects, as well as other effects such as the lunar cycle 
and their effects on the photometry, we instead choose a simpler 
prescription. Each event is monitored with continuous photometry 
at a reduced cadence of 30 minutes. These observations are per- 
formed by telescopes with 1.3 m effective diameter observing in 
the /-band. For each exposure of 120 s, the seeing is chosen from a 
lognormal distribution with mean 1.2 arcsec and standard deviation 
0.25 arcsec, and a background flux distributed as 



F = 8500LN(1.5,0.4) photons arcsec" 



(15) 



which is integrated over a seeing disk, and where LN(^i, a) is a 
lognormal distribution with mean /j, and standard deviation a. New 
values of seeing and background flux are chosen for each obser- 
vation. A lower limit on the photometric accuracy is imposed by 
adding a Gaussian noise component, with dispersion 0.3 percent, 
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Figure 3. An example lightcurve from the simulations that required cov- 
erage to be extended to cover a feature far from the lightcurve peak. The 
upper panel shows the lightcurve (/-band magnitude) and the lower panel 
shows the Paczyriski residual (/-band residuals from the single lens fit). Red 
points show the simulated data points with error bars, and black, green and 
blue lines are the best-fitting Paczyriski model, the best-fitting static binary 
model and the true orbital motion model (largely hidden below the green 
static model curve) respectively. Only one data point in 24 is shown for 
clarity. The lightcurve shown is for that of an event by a stellar binary lens 
with q = 0.22, d w 8.6 and t E = 14.9 d. Usually, only data points that 
cover the inner 54e are used, apart from some data points used to constrain 
the baseline magnitude (cf. the lightcurve for times t > d), however ad- 
ditional data points are used to fully cover the additional lightcurve feature 
down to baseline (cf. the lightcurve for t < d). Further details for the 
event can be found in Tables 4 and 5 in the online supplementary material. 



to the photon counts, which are calculated by assuming 10 pho- 
tons m -2 s _1 reach the observer from a / = 22 source. 

To ensure that all the features of a lightcurve are covered, and 
that there is a good balance between the baseline, peak and features 
of the lightcurve when fitting (see the next section), the lightcurve 
is monitored continuously over the times —5fE ^ t— to < 5*e, and 
over 10.5tE ^ \t — to| < 9.5£e to sample the baseline. To ensure 
that all features are covered, if the magnification of the source rises 
above A ^ Athrcsh = 1.0062, the coverage is extended so as 
to be continuous within one Einstein time-scale of the feature and 
continuous between the feature and t — to. Figure [3] shows an 
example of a lightcurve where coverage had to be extended. 



4 MEASURING ORBITAL MOTION 

Ultimately we are interested in finding the fraction of binary mi- 
crolensing events which show signs of orbital motion. To do this 
we must classify the events we simulate into those binary events 
that do show orbital motion, those that do not, and events that do 
not show binary signatures. To do this we fit each event first with a 
single lens model, and then those events which are poorly fit with 
this model, we fit with a static binary lens model. To evaluate the 
effectiveness of each stage of the fitting process, in addition to the 
fitting of the lightcurves simulated with orbiting binary lenses, we 
must also fit lightcurves simulated with point-mass lenses and static 
binary lenses. 



4.1 Lightcurve modelling 

The single lens model has five parameters: the time of closest ap- 
proach to, the event time-scale t E , the impact parameter , the 
baseline magnitude if and the blending fraction ff . We performed 



a x minimization using the MINUIT routine from CERNLIB (James 



|& Roos||l975| >, with all parameters free; all parameters were un- 
constrained, except for ff, which was constrained to be within 
0-0 < ff < 1.2. For each event, seven single lens fits were per- 
formed, with different initial blending fractions, ff — 0.05, 0.2, 
0.4, 0.6, 0.8, 1.0 and 1.2. For each fit, the initial guesses for each 
parameter were: tf — 0, the time-scale was the true time-scale, 
the baseline magnitude was taken to be the magnitude of the first 
data point on the lightcurve, and the impact parameter was chosen 
such that at t — to the magnitude of the event would be that of the 
brightest data point. This prescription works well for events which 
are well modelled by a single lens model, but not so for events 
with strong binary features, or events which are heavily blended 
and barely rise above baseline. It is therefore useful to eliminate 
events falling into the later category before performing the fitting, 
such that the only events that the single-lens model fails to fit are 
ones that show genuine signs of lens binarity. This cut will be de- 
scribed in the next subsection. 

To fit the binary lens lightcurves, we found it necessary to split 
the events into caustic-crossing events and non-caustic-crossing 
events, and to fit each category using a different parametriza- 
tion. The non-caustic crossing events we fitted with a standard 
parametrization, with a reference frame centred on the primary 
lens. The parameters are: the time of closest approach to the lens 
primary to, the event time-scale t%, the impact parameter between 
the lens primary and the source Uo, the angle of the source trajec- 
tory to the binary axis 8q, the logarithm of the projected binary 
separation log d s , the logarithm of the normalized secondary mass 
logm^, the baseline magnitude if and the blending fraction ff. 
For brevity we introduce the vector notation 



(t s o,t%,ul6l,\o g d s ,logmi,/f, ff) , (16) 



to represent the parameter set of the standard binary parametriza- 
tion. 

For the number of lightcurves necessary to obtain a good sta- 
tistical sample, a full search of the full binary lens parameter space 
is not computationally feasible, so we perform just one minimiza- 
tion per lightcurve. We must therefore pay special attention to the 
choice of initial guesses we use, firstly so as to maximize the chance 
of finding a good minimum, and secondly so as to treat the fitting 
of the static binary events comparably to the orbiting binary events. 
The static binary simulations are drawn from the same distributions 
as the orbiting binary simulations, the only difference being that the 
lens is frozen in the state it would be in at t = to. 

As we have simulated the microlensing events, we already 
have a perfect knowledge of the systems, and we can use this 
knowledge to obtain a good set of initial guesses. We note that at a 
given time, the state of an orbiting binary lens can be described by 
a static binary model. We can therefore describe our lens at time t 
using the time dependent parameter set 



p(t) = (to,t E ,u ,Oo(t),d(t),q,I b , f s 



(17) 



where we have used the same definitions and centre of mass ref- 
erence frame as in the previous section. Note that only two of the 
parameters are time dependent, and so we can use the true values of 
the constant parameters as initial guesses, having applied the appro- 
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priate coordinate transformations [J However, we are still left with 
the problem of choosing the guesses of 8q and d s . We could choose 
0o (to) and d(to), but this would bias the fitting success probability 
unfairly towards static binary events: the initial guess would be the 
actual model used to simulate the data. Instead, we choose to use 
d(tf ) and 8o(t{), where U is the time of a feature in the lightcurve. 
We define a feature simply as any maximum in the lightcurve, or 
a maximum or minimum in the Paczynski residual (the residual of 
the true lightcurve with respect to the best fitting single lens model) 
with \ I — Jpac| > 0.1, where / is the J-band magnitude of the true 
model, and Jp ac the /-band magnitude of the best fitting Paczynski 
model. As there are in general more than one feature, we choose the 
feature that gives the best x 2 (p(if))- If the initial guesses for fits to 
static binary lightcurves are chosen in the same way, as if the binary 
were orbiting, then the initial guesses for static lenses should be 
worse than for orbiting lenses, as at the time of the chosen feature, 
the true orbiting lens magnification will exactly match the magnifi- 
cation of the initial guess static model. In reality, for ti ~ to there 
will likely be a bias in favour of static lenses and tf 56 to there 
will be a bias in favour orbiting lenses, but we do not believe this 
will affect results significantly. To fit the events, we again use the 
MINUIT minimizer, allowing all parameters to vary. All parameters 
are unconstrained, except for ff , which is constrained to the range 
< ff < 1.2. 

While this method was suitable for events which showed 
smooth binary features, it is not always suitable for those events 
which exhibit caustic crossings. For these events, in addition to fit- 
ting with the standard parametrization, we also used the alterna- 
tive parametrization of Cassan (2008). This replaces the parame- 

' ' ' ' Q Q G Q 

ters specifying the source trajectory (to, £e, Uq, oq), with param- 
eters that better reflect the sharp caustic crossing features of the 
lightcurve (i^ n , i^? x , s ex) the times of a caustic entry and exit 
and the positions of the entry and exit on the caustic respectively, 
where sj? n and s^ x , are defined to be the chord length along the 
caustic, normalized such that ^ sj? n ) s cx < 2. Full details of 
the parametrization can be found in the Cassan ( 2008 1 paper. The 
parameter set we use for caustic-crossing events is therefore 

P° = s2»s2o lo S rfC , lo S q°,Ib,f? J , (18) 

where the parameter logm^ has been replaced by logg c as a 
matter of preference; the two parameters are related by mf = 

g c /(i + g c ). 

The accurate calculation of the sj? n and s^ x parameters is 
quite computationally expensive, compared to the calculation of 
a lightcurve, and needs to be repeated each time d or q changes. 
Also, despite the improved parametrization, the \ 2 surface is still 
very complicated, especially in the s£,~ s ex plane, containing many 
local minima. For these reasons we pursue a three stage mini- 
mization process. We begin by conducting a grid search over the 
entire s^-s^ plane, with 128 x 128 points spaced evenly in 
s en> s cx> an d with all other parameters, including the caustic cross- 
ing times, fixed at their true values, except for logd c . logd c is 
fixed at a random value is chosen in from the range A log d c — 
1.5(logd(t e x) — logd(ten)) or Alogd = 0.015, whichever is 
greater, centred on the midpoint of log d between the caustic cross- 
ings, and where d(t en ) and d(i cx ) are the projected separations at 
the caustic entry and exit times respectively. The range of A log d c 



3 In the reference frame of , to and Uo would also be time dependent as 
the origin (the primary mass) is not fixed. 
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Figure 4. Example lightcurve and caustic map of an event where a non- 
adjacent caustic entry-exit pair was chosen for fitting with the |Cassan| 
1 2008! parametrization. The lightcurve is shown in the left panel, where 
red points show the simulated data, the blue line is the true model and the 
green line is the static binary model. The right panel shows a map of the 
caustic of the static binary model, plotted in red, and the source trajectory, 
plotted in black. The numbers indicate the order of the caustic crossings. 
The static model has been adjusted by hand to aid clarity. Further details 
for the event can be found in Tables 4 and 5 in the online supplementary 
material. 



is truncated, if necessary, to ensure that it only covers the caus- 
tic topologies at the time of the crossings. For the static lenses, 
log d c is chosen from a uniform distribution with the same range 
as if the lens were orbiting. The grid search is then refined by per- 
forming a second 128 x 128 grid search over a box of side length 
1/32 about the grid point with the lowest \ 2 - Six 2 x 128 x 128 
grid searches are performed with different random values of log d c . 
In cases where there are multiple caustic crossings, different pairs 
of caustic crossings are used to define (t^ n , t^, s^ n , s^) for each 
grid search. Figure [4] shows an example lightcurve where the first 
caustic exit defines (t^ n , sj n ) and the second caustic entry defines 
(t c s c ) 

Vex ) "ex/ • 

The second stage of the fitting simply polishes the result of 
the grid search by performing a minimization over s^ n and sj? x , 
with all other parameters fixed, using MINUIT. In the final stage 
of the fitting, all parameters, except for t^ n and tj^ are allowed to 
vary in a further MINUIT minimization. Again, all parameters were 
unconstrained, except for which was constrained to the range 
< /° < 1.2. We found that, at all stages of the minimization for 
caustic crossing events, the minimization performed better when 
the first and last data points inside the caustic crossing were not 
considered in the fit. This is because, with the high cadence obser- 
vations that we simulate, the point source is typically very close 
to the inside of the fold caustic, and hence is magnified by many 
orders of magnitude. This leads to unrealistic photometry in two 
ways: firstly, in a real detector, saturation would become a prob- 
lem, and secondly, a real, finite, source would not be magnified in 
such an extreme way. 



4.2 Classification of events 

With the modelling procedures in place, we now describe the clas- 
sification of the events. The classification is performed by a series 
of cuts based on the % 2 results of the fitting described in the last 
subsection. The first cut, the variability cut, removes events which 
do not show significant variability from the analysis. This is done, 
without fitting, by comparing the \ 2 values of the simulated data 
relative to the true model, Xom» an d relative to a constant lightcurve 
with no variability at the true baseline magnitude, xi- We exclude 
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events that do not satisfy 



Table 1. Summary of the results for planetary lenses. 



«obs 



xl 



2 

XOM 



> 0.3, 



^obs 



(19) 



where n D b s is the number of observations. 

The second cut is used to classify events into single lens-like 
events, and binary lens-like events, or events that do not and do 
exhibit binary lens features in their lightcurves. Using the results of 
the single lens modelling, xpao me X 2 °f me simulated data with 
respect to the single lens model, we define events that satisfy 



AXPac 



2 

XPa 



Xom > 200 



(20) 



to be binary events, and those that do not to be single events. Bi- 
nary events can then be split into caustic crossing binary events and 
smoothly varying events, or caustic crossing and smooth events re- 
spectively. We define a caustic crossing event as one where at least 
one data point is measured when the source is inside a caustic]^] 

The final cut is based on the result of lightcurve fitting with 
binary models. Events that satisfy 



Ax! 



Xl - Xom > 200 



(21) 



Orbit 


static 


circular 




Single 


48511 


49226 




Binary 


1364 


1366 




Caustic 


410 


449 




Caustic static 


397 


414 




Caustic orbital motion 


7 


35 




Smooth 


954 


917 




Smooth static 


931 


883 




Smooth orbital motion 


23 


34 




Table 2. Summary of the results for stellar binary lenses. 


Orbit 


static 


circular 


eccentric 


Single 


4151 


4046 


4153 


Binary 


1413 


1424 


1385 


Caustic 


641 


635 


613 


Caustic static 


608 


538 


550 


Caustic orbital motion 


25 


86 


61 


Smooth 


772 


789 


772 


Smooth static 


764 


743 


729 


Smooth orbital motion 


g 


46 


43 



are classified as events that exhibit orbital motion (orbital motion 
events) and those that do not are classified as static events, where 
Xs is taken to be the \ 2 of the best fitting static binary model. 
For smooth events this is the \ 2 of the best fitting standard binary 
model, and for caustic crossing events it is the \ 2 of the better 
fitting of the Cassan ( 2008 1 caustic crossing model or the standard 
binary model. In the case of the caustic crossing fits, the data points 
removed from the lightcurve do not contribute to Xom- 

With these classifications in place, we can now define the 
binary detection efficiency and the orbital motion detection effi- 
ciency. The binary detection efficiency is the fraction of detectable 
microlensing events that show binary signatures 



CBS 



Nbs 



(22) 



where N m i is the number of events satisfying Axb/ n obs > 0.3, 
and Nbs is the number of events satisfying Axp ac > 200. The 
orbital motion detection efficiency is the fraction of binary events 
that show orbital motion signatures 



£OM 



A^om 



(23) 



5 RESULTS 

5.1 What fraction of events show orbital motion? 

We begin by presenting and analyzing the results of the simulations 
as a whole, calculating the fraction of microlensing events in which 
we expect to see orbital motion events. Tables [T] and [2] summarize 
the results of the cuts described in the previous section, for plane- 
tary and stellar binary events respectively. It should be noted that 
in a small number of caustic crossing events, the fitting procedure 
failed, and so these events have been excluded from the analysis of 
the orbital motion detection efficiency, but not of the binary detec- 
tion efficiency. These events are included in the Binary and Caustic 
rows of Tables [T] and [2] but not in the others. Figure [5] shows some 
lightcurves which were slightly below the threshold for each cut. 

Table [3] shows the binary detection efficiency and orbital mo- 
tion detection efficiency for both planetary and stellar binary lenses. 
It should be noted that the binary detection efficiency will be larger 
than for microlensing events with finite sources, as the effect of 
the finite source will be to smooth out sharper lightcurve features, 

4 The removal of data points in the fitting process does not affect the clas- 
sification. 



where ./Vom is the number of events satisfying A%s > 200. 

To be confident of our results we must quantify the effec- 
tiveness of the modelling prescriptions we use. We can do this by 
measuring the rate of false positives in our samples. To measure 
these rates we simulate both single lens events and static binary 
lens events, drawn from the same distributions as the orbiting lens 
events. These events then go through the same fitting procedure as 
the orbiting lens events and are subject to the same cuts. The binary 
lens false positive rate £Bg gl ° is therefore the fraction of detectable 
single lens microlensing events that survive the Axlac > 200 cut, 
and the orbital motion false positive rate £om c ls me fraction of 
static binary lens events that survive the Axl > 200 cut. 



Table 3. Binary and orbital motion detection efficiencies. 



Orbit 




circular 


eccentric 


q < 0.01 


e BS 


0.0772 ± 0.0014 




q < 0.01 Caustic 


e OM 


0.061 ± 0.010 




q < 0.01 Smooth 


e OM 


0.0130 ±0.0055 




q < 0.01 All 


e OM 


0.029 ± 0.005 




q > 0.01 


CBS 


0.260 ±0.004 


0.251 ± 0.004 


q > 0.01 Caustic 


e OM 


0.098 ±0.011 


0.060 ±0.010 


q > 0.01 Smooth 


e OM 


0.048 ± 0.006 


0.045 ± 0.006 


q ^ 0.01 All 


£ OM 


0.070 ± 0.006 


0.052 ±0.006 
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Figure 5. Example lightcurves of three events that narrowly failed one 
of the classification cuts. From top to bottom, the lightcurves failed the 

A Xb _ Xb~XQM > g g ^^,2 _ ^,2 _ ^,2 



2 

Pac 



Xpa 



Ax 2 



Xqm > 200 ' ^ 

s = x| — Xom "^00 cuts, respectively. The upper panel of each 
subplot shows the lightcurve, and the lower panel the residual with respect 
to the appropriate model for the cut, i.e. the constant baseline 'model', the 
best fitting Paczyhski model, and the best fitting static binary model in the 
top, middle and lower subplots, respectively. Colour coding is the same as 
in Figure]!] and the cyan line in the top subplot shows the constant baseline 
'model'. Further details for the events can be found in Tables 4 and 5 in the 
online supplementary material. 



and usually reduce the amplitude of deviations from the single 
lens model. This means that ees for planetary lenses is likely a 
significant overestimate, however, for stellar binary lenses the re- 
sult is likely to be more realistic as binary lightcurve features tend 
to be stronger and have longer durations. The detection efficien- 
cies presented have been corrected for systematic false positives 
from each fitting stage by subtracting the measured false positive 



positive rates of £om 1c 



static 



for orbiting lenses. From a simulation of 10 single lenses with no 
false positives we measured e^g 810 = 01q 7x10 , where the error 
quoted is a statistical la confidence limit calculated using Wilson's 
score method (Wilson 1927, Newcombe 1998b). To calculate the 
errors on the corrected detection efficiencies shown in the table, 
and on those we present in the next subsection, we use Wilson's 
score method adapted for the difference of two proportions JNew-| 
|combe|1998a| method 10). For planetary events we measured false 
0.0241+°;°°^ for smooth events and 
o 0039 f° r caus ti c crossing events. For stellar 
binary events we measured £om ic = 0.0104t ; u0 22 f° r smooth 
events and £ " ic = 0.0395±o'oo50 for cau stic crossing events. 
The overall orbital motion detection efficiencies were calculated 
as a weighted average of the detection efficiencies for smooth and 
caustic crossing events, once corrected for false positives. 

While in many cases we may not be able to say that a 
lightcurve in our simulations definitively shows orbital motion sig- 
natures, due to relatively high rates of false positive detections, 
there is a clear excess of detections in the circular and eccentric 
orbit simulations relative to the static ones, though the detection of 
this excess is only marginal in smooth planetary events. Interest- 
ingly, there appears to be a discrepancy in the orbital motion detec- 
tion efficiencies for stellar binary caustic crossing events. The same 
static orbit simulation results were used to calculate the corrected 
orbital motion efficiencies for both circular and eccentric orbits, 
which means that the measurements are not independent. Also, the 
eccentricity of the orbits allows the projected separation to take 
a wider range of values than the circular orbits, which means the 
false positive rate measured with the same distribution for circu- 
lar orbits is likely an overestimate for eccentric orbits; for caustic 
crossing events the majority of false positives are caused by events 
with resonant caustic topology (see Figure[T9]later in this section). 
We therefore believe the discrepancy to be caused largely due to a 
combination of a relatively large statistical fluctuation in the num- 
ber of eccentric orbit events that do show orbital motion, and an 
overestimate of the false positive rate for eccentric orbits. 



5.2 What affects the detectability of orbital motion? 

We now investigate the effects that various system parameters have 
on the detectability of orbital motion. We look at the dependence 
of the orbital motion detectability on both the standard microlens- 
ing parameters and the physical orbital parameters, and compare 
them where appropriate. While we conducted two sets of simula- 
tions, one with circular orbits and one with eccentric orbits, we only 
present the results for those with circular orbits here, as both sets 
are in good agreement. 

We begin by looking at the dependence on the impact pa- 
rameter, the sole parameter that determines the maximum mag- 
nification of a sin gle lens rnicro lensing event j4 max = (uq + 
2)/(u -\Ao + 4 ) i Paczynski|l986 i. For all binary lenses, except 
wide stellar binaries, the central caustic is located near to the cen- 
tre of mass, and so uq determines whether or not the source will 
encounter this caustic. Figure [6] plots the orbital motion detection 
efficiency as a fraction of caustic crossing or smooth binary events 
(top panels), and the total number of orbital motion detections (bot- 
tom panels), against the impact parameter for both planets (left pan- 
els) and binary stars (right panels). In the plots, red lines represent 
data for caustic crossing events and blue lines for smooth events. In 
the top panel the orbital motion detection efficiency has been cor- 
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Figure 6. Plot of the orbital motion detection efficiency, corrected for sys- 
tematic false positives (top panels), and the absolute number of orbital mo- 
tion detections in the simulations (lower panels), against the impact param- 
eter |no|. Results are shown for lenses with planetary mass ratios (left) and 
binary star mass ratios (right). Red lines with filled squares show the results 
for caustic crossing events and blue lines with filled circles show the results 
for smooth events. In the upper panels a line marks zero orbital motion de- 
tection efficiency. All events had circular orbits, and in the lower panels 
results are shown for events where the lens components were in orbit (solid 
lines, filled points) and where they were held static for the calculation of the 
false positive rate (dashed lines, open points). Events have been binned into 
bins of equal width, and points plotted at the centre of the bin. Note that in 
the lower panels the scales are different. 



rected for systematic false detections as described in the previous 
subsection, whereas the bottom panel shows the number of detec- 
tions for both orbiting (solid lines, filled points) and static lenses 
(dashed lines, open points). Note that the orbital motion detection 
efficiency can be negative, due to statistical fluctuations, but if it is, 
the measurement should be consistent with zero. The events have 
been binned into bins of constant width, on the scale that they are 
plotted. It should also be noted, that the number of planetary events 
simulated was a factor of 9 larger than the number of stellar binary 
events. 

The plots of orbital motion detection efficiency (from here on 
detection efficiency) against |«o| for caustic crossing events show 
much the same trends for both planetary and stellar binary lenses, 
with significant detection efficiencies for high-magnification (low 
| Mo |) events only, with no caustic crossing planetary detections for 
|«o| > 0.6, and only a few for stellar binaries. This is due to the lo- 
cation of central and resonant caustics close to the center of mass, 
which the source can only cross in events with small |ito|. Con- 
sequently, for the events with larger uo, the source can only cross 
weaker secondary caustics, which in the case of wide binaries will 
typically move slowly, and in the case of close binaries are typi- 
cally very small and are rarely crossed. The secondary caustics of 
close stellar binaries are significantly larger and stronger than those 
of planetary lenses, and so the chances of the source crossing them 
is higher, and the caustic has a longer time in which to change due 
to orbital motion as the source crosses it, leading to the small posi- 
tive efficiency for \uo\ > 0.6. For smooth events, the planetary and 
stellar binary lenses show weak but opposing trends, with the effi- 
ciency increasing slightly as \uo \ increases for planetary events and 
decreasing slightly as |ttn| increases for stellar binary events, indi- 
cating that the impact parameter only plays a small role in orbital 
motion detectability for smooth lightcurves. Note however, that for 
both smooth and caustic crossing events the number of orbital mo- 
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Figure 7. As Figure|6] but plotted against the event time-scale tg. 



tion detections, as opposed to the detection efficiency, is a strong 
function of |uo|, peaking at small values, due to the dependence of 
the binary detection efficiency on the impact parameter. 

Figure [7] plots the detection efficiency against the event time- 
scale £e- All classes of binary event (planetary or binary, smooth 
or caustic crossing) show a strong detection efficiency dependence 
on the event time-scale. The reason for this dependence is simply 
because a longer time-scale allows the lens to complete a larger 
fraction of its orbit, and hence cause a larger change in the magni- 
fication pattern, during the time in which the source probes regions 
of the magnification pattern that deviate from that of a single lens. 
In the case of planetary lenses, it seems that a time-scale of greater 
than ~ 10 days is necessary for caustic crossing events and slightly 
longer for smooth events. Caustic crossing events show larger de- 
tection efficiency than smooth events, even at shorter time-scales. 
This is likely due to the high accuracy with which caustic cross- 
ing times, and the lightcurve shape around caustic crossings can be 
measured. In the case of OGLE-2006-BLG-109 this has allowed 
the orbital motion of the lens to be measured from data covering 
just ~ 0.2 percent of the orbit ( |Gaudi et al.||2008[ [Bennett et al.] 
|2010| >. Smooth events in contrast require a much larger fraction of 
the orbit to cause significantly detectable changes in the lightcurve, 
and hence require a longer time-scale to achieve the same detection 
efficiency. However, typically it is possible for smooth features to 
cover a much larger fraction of the lightcurve than caustic crossing 
features, lessening the effect of this discrepancy. 

For stellar binary lenses, orbital motion features can be can be 
detected effectively over almost the entire range of time-scales that 
we simulated, though with a low efficiency for time-scales below 
~ 40 days for smooth events and ~ 10 days for caustic crossing 
events. For events with time-scales over ~ 100 days, the detection 
efficiency reaches ~ 20 percent for smooth events and ~ 40 per- 
cent for caustic crossing events. The detection efficiencies are simi- 
lar for planetary events. The majority of planetary and binary events 
showing orbital motion have time-scales of around ~ 10-40 days, 
with few events at larger £e due to the steep t^ 3 distribution at large 
time-scales (Mao & Paczynski 19961. However, the strong depen- 
dence of eoM on time-scale means that the slope of the high t-E tail 
of the distribution of orbital motion events is much shallower than 
t -3 

The plots of detection efficiency against projected separation 
do (Figure[8]l and semimajor axis a (Figure[9) tell largely the same 
story. The detection efficiency in stellar binaries has a significant 
inverse dependence on both do and a, as would be expected from 
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Figure 8. As Figure[6] but plotted against do the lens separation at time to . 



-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 
log q logq 

Figure 10. As Figure[6] but plotted against the mass ratio q. 
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Figure 9. As Figure[6] but plotted against the semimajor axis a. 

the dependence of orbital velocity of semimajor axis. However, the 
behaviour for planetary lenses is less intuitive: for caustic cross- 
ing events, there is a significant peak in the detection efficiency at 
a ~ 4 AU, and a peak/shoulder at do ~ 2. There is a second peak 
in eoM with do. The two peaks occur at values of do where the 
boundaries between caustic topologies occur for the highest mass 
ratio planets. It is at these boundaries that, for a small change in 
projected separation d(logd), the largest changes in the caustics 
occur. The peak in com against a at a ~ 4 AU for caustic crossing 
planetary events is accompanied by a hint of a peak at small values 
of a. The peak at a ~ 4 AU can be explained by considering the 
typical scale of the Einstein ring, and by considering the trend of 
eoM with the event time-scale. The typical size of the Einstein ring 
for a microlensing event is 2-3 AU, but as seen in Figure [7] orbital 
motion effects typically occur in events with larger time-scales. As 
the time-scale is correlated with the Einstein ring size, and caustic 
crossing events typically occur in systems with do ~ 1, the peak or- 
bital motion detection efficiency occurs at a semimajor axis slightly 
above the typical Einstein ring size, at a ~ 4 AU. The increase in 
orbital velocity as a decreases likely causes the second weaker peak 
in eoM at smaller a. Little can be said about the trend of eoM with 
a for smooth planetary events as small numbers of events, and the 
distribution of Einstein radius sizes serves to smear out any obvi- 
ous trends. However, when plotted against do, eoM does increase 
towards smaller values of do, as would be expected from orbital 
velocity considerations. 



Returning to the caustic crossing stellar binary events, com 
flattens off as a increases to ~ 4 AU, before dropping to zero. 
This flattening likely has the same cause as the peak for planetary 
caustic crossing events. We see the more intuitive inverse trend in 
stellar binaries because of the stronger and larger magnification pat- 
tern features that they exhibit, and the larger range of d over which 
the caustics have a significant size. This results in a distribution of 
events over a and do which is broader and somewhat less peaked 
than for planetary events (see the lower panels of Figures[8]and|9]l. 
This allows the inverse relationship between orbital velocity and 
semimajor axis to have a greater influence on the trend in the or- 
bital motion detection efficiency. We note that the reason we see 
such a complicated relationship between com and a and do, but 
not for example between com and ts,, is that the orbital separation 
affects the orbital velocity in a relatively simple way and caustic 
size and strength in a complicated way, whereas ts only affects, or 
more accurately is the result of, a fairly simple dependence on a 
single factor in the detection of orbital motion, the source speed. 

Figure [10] plots the detection efficiency against the mass ra- 
tio q. Treating both planetary and stellar binary lenses together, 
there is a trend of increasing detection efficiency with increasing q, 
for both smooth and caustic crossing events. However, for caustic 
crossing events, this increase is very shallow, with a factor of < 3 
increase over three decades in q, from logg « —3 to logg = 0. 
For smooth events, there is a stronger trend, with the detection ef- 
ficiency being effectively zero for log g < —3.5, while rising from 
~ 1 percent to ~ 10 percent over the range —3.5 < logg < 0. 
These shallow dependencies are somewhat unexpected in relation 
to the somewhat stronger g 0,5 dependence of the binary detection 
efficiency, which derives directly from the dependence of caustic 
size on g (Han 2006). However, the orbital detection efficiency ef- 
fectively divides through by this dependence (unlike the curves of 
the number of orbital motion detections, which show a strong de- 
pendence on g), to leave a very shallow orbital motion detection 
efficiency curves. The other effect that g has on the lightcurve fea- 
tures is to make them stronger as q increases. In caustic crossing 
events the caustic features are usually strong, independent of the 
value of g, and hence the caustic crossing events curve is shallower 
than the curve for smooth events, for which the dependence of the 
feature strength on q is much more important. 

Figure [TJJ shows the detection efficiency plotted against the 
primary lens mass. The dependence is as expected for both mass 
ratio regimes and for both types of binary event, increasing as the 
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Figure 13. As Figure|6] but plotted against the orbital period T. 
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Figure 14. As Figure|6] but plotted against the baseline magnitude 7b- 



mass of the primary increases. The trend is strongest in smooth, 
stellar binary events. 

Figure [T2] plots the detection efficiency against the lens dis- 
tance. In all cases, a trend of increasing detection efficiency with 
decreasing lens distance is seen, though caustic crossing events 
suffer from small number statistics at low values of D\/D s . Note 
however, that the frequency distribution (plotted in the lower pan- 
els of Figure [T"2) of orbital motion events, once false positives have 
been approximately accounted for, is different, being peaked at 
Di/Ds ~ 0.7. 

Figure [T3] shows the detection efficiency plotted against the 
orbital period. Both types of stellar binary event show a significant 
inverse trend. Planetary caustic crossing events show a peak, and 
stellar caustic crossing events a flattening, at large periods. These 
features correspond directly to similar features in the curves of com 
with a and will have the same cause. 

Figures [14] and [T5] plot the detection efficiency against the 
baseline magnitude lb and blending fraction / s respectively. For 
our purposes, the primary effect of both parameters is to affect the 
accuracy with which microlensing variations can be measured in 
the lightcurve. For a fixed observing setup, the baseline magni- 
tude determines the photometric accuracy, which should lead to a 
trend of increasing detection efficiency with decreasing magnitude. 
This is seen to a certain extent in all cases, but brighter events may 
suffer significantly from blending, due to faint source stars falling 
entirely within the large point spread function of a much brighter 



star. Blending determines the relative strength of features in the 
lightcurve, and as such has a much more significant effect on the de- 
tection of smooth binary features, which have a continuous range of 
shapes and sizes, compared to the effect on caustic crossings which 
are typically sharp and very strong, at least when finite sources are 
not considered. It is no surprise, therefore, that smooth stellar bi- 
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Figure 15. As Figure [6] but plotted against the fraction of baseline flux 
associated with the source f s . 
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Figure 16. As Figure[6] but plotted against the orbital inclination i 



nary events show a significant increase in orbital motion detection 
efficiency with blending fraction. This is less obvious in planetary 
lenses, likely because the smooth lightcurve features of planetary 
lenses are often very weak and difficult to detect even without the 
hindrance of the blending, and would not permit the measurement 
of higher order effects for any value of blending fraction. It is more 
surprising, perhaps, that caustic crossing events show a significant 
dependence on blending, as in the simulations all caustic crossing 
events were detected as binaries, regardless of blending. This im- 
plies that, at least in some orbital motion detections in caustic cross- 
ing events, the additional smooth features in the lightcurve, such as 
peaks and shoulders due to cusp approaches outside the caustic, and 
features due to fold caustic approaches within the caustic, play an 
important role in the detection of orbital motion (e.g. lightcurves a 
and e in Figure 2j]in the next subsection). 

Figure [16 plots the detection efficiency against inclination. 
There is little evidence for any significant dependence on inclina- 
tion for all caustic crossing events, and for smooth planetary events. 
There is however a stronger trend for smooth stellar binary events, 
the detection efficiency decreasing as the inclination increases. This 
would be expected in systems where cl/ve < d c , the boundary be- 
tween close and resonant caustic topologies, where a reduction in 
the projected separation due to inclination would reduce the size 
of the caustics and reduce the detectability of both binary features 
and orbital motion signatures. Unfortunately, due to the similar ef- 
fects of inclination and eccentricity on the projected orbit, the data 
from the eccentric orbit simulations did not show any dependence 
of eoM with eccentricity. This however implies that the effects of 
eccentricity on the orbital motion detection efficiency are not likely 
to be significantly stronger than those of inclination. 

It is important not just to consider the system parameters in 
isolation, but also their combined effects on the orbital motion de- 
tection efficiency. For example, Dominik ( 1998 1 introduced two di- 
mensionless ratios to describe the magnitude of orbital motion ef- 
fects on a binary lens: 



T' 



the ratio of time-scales, and 



R v — 



v t 



(24) 



(25) 



the ratio of velocities, where v c i IC = a/2nT is the circular velocity 
of the orbit. These ratios attempt to encapsulate the most important 
factors that determine if an event will show orbital motion features. 
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Figure 17. As Figure [6] but plotted against the ratio of microlensing to 
orbital time-scales Rt =t^ J /T 
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Figure 18. As Figure [6] but plotted against the ratio of orbital and source 
velocities R v = i) c i rc / vt 



Figures^]and^]plot the detection efficiency against Rt and R v 
respectively. Both ratios prove to be good descriptors of the orbital 
motion detection efficiency, with com showing strong increasing 
trends as Rt and R v increase, across all mass ratios and lightcurve 
types, though with a lower significance in planetary events. It would 
even seem that, in the case of smooth events, there exists a thresh- 
old value of the ratios, below which the orbital motion detection 
efficiency is negligible. For the ratio of time-scales, the threshold is 
log Rt ~ — 2 for both planetary and stellar binary lenses, while for 
the ratio of velocities the value appears to be more dependent on the 
mass ratio, taking values of logf? v ~ —2.5 for planetary lenses, 
and log R v « —2.75 for stellar binary lenses. There may be sim- 
ilar thresholds for caustic crossing events, but at smaller values of 
Rt and 7? v . 



5.3 Are there two classes of orbital motion event? 



|Gaudi| ^2009} has suggested that orbital motion can affect the 
lightcurves of microlensing events in two ways. In the first sce- 
nario, the orbital motion effects are dominated by rotation in the 
lens, as the orientation of binary axis changes during the time be- 
tween two widely separated lightcurve features. The second type of 
effect is due to changes in the projected separation over the course 
of a single lightcurve feature such as a resonant caustic crossing. In 



14 M. T. Penny, S. Mao and E. Kerins 



this subsection we will describe the typical features of each type of 
event before investigating to what extent orbital motion events can 
be classified in such a way. 

Gaudi (20091 describes the separational class of event as typ- 
ically occurring in archetypal binary microlenses with resonant 
caustic crossings. If the binary's orbit is inclined, the projected 
separation of the lenses changes, causing a stretching or compres- 
sion of the resonant caustic. If the projected separation is close to 
a boundary between caustic topologies, d ~ d c , or d ~ d w , the 
changes in the caustic structure can be very rapid. If the microlens- 
ing event occurs while the changes are happening, and the source 
crosses, or passes close to, the caustics, there is a very good chance 
of detecting the orbital motion. As a whole though, the changes 
in caustic structure during the caustic crossing time-scale will be 
fairly small, e.g. the difference in caustic crossing time between the 
static lens and the orbiting lens may be of order minutes to hours 
(cf. the orbital period of several years). It is only the extremely good 
accuracy with which caustic crossings can be measured and timed 
that facilitates the high orbital motion detection probability. These 
changes to the caustic shape will often be more significant than the 
changes in orientation of the caustic due to rotation, and so we class 
them as separational orbital motion effects. 

|Gaudi|p009| > described the rotational class of event as occur- 
ring when a source encounters two disjoint caustics of a typically 
close topology lens. In the time between the two caustic encounters, 
which are separated by a time At ~ £e, the lens components have 
time to rotate and show detectable signatures of orbital motion. We 
extend the class by considering the important effect to be the long 
baseline over which binary lensing features can be detected. If bi- 
nary lens features are detectable across a significant fraction of the 
lightcurve then a significant amount of rotation can occur in the 
lens while the features are detectable. Up to now, our discussion 
has focused mainly on caustic features, whether the source crosses 
them or not, but, in stellar binary lenses especially, the magnifica- 
tion pattern of the lens can differ significantly, if subtly, from the 
single lens form over large parts of the pattern, and well away from 
caustics. For example, in close binary lenses, there is a region of ex- 
cess magnification that can stretch the entire distance between the 
facing cusps of the central and secondary caustics. In stellar binary 
lenses, this can extend for distances larger than an Einstein radius. 
In planetary lenses the magnification excesses are weaker, but there 
tends to be a large region of demagnification between the two plan- 
etary caustics. If lenses with such features rotate rapidly, then the 
source may encounter them in such a way that a static lens inter- 
pretation of the lightcurve features is not possible, and lens rotation 
must be invoked. 

We begin by looking for evidence of two classes of event in 
the locations of the orbital motion events in the do-q plane. Fig- 
ure[l9]plots q against do for all binary events; events which do not 
show orbital motion signatures are plotted with small, open points 
with light colours, whereas those that do are plotted with large, 
filled points with darker colours. Caustic crossing events are plot- 
ted with red squares, and smooth events with blue circles. Upper 
panels show stellar binary lenses and lower panels show planetary 
lenses, while the left panels show orbiting lenses and the right panel 
show static lenses. The black lines show the boundaries between 
the caustic topologies (equations [7] and [5}. It is immediately clear 
that caustic crossing and smooth orbital motion events reside in 
different regions of the do-q plane, with virtually all events within 
the intermediate topology regime being caustic crossing. Almost 
all smooth orbital motion events are located in the close topology 
region. This broadly reflects the underlying pattern for all binary 
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Figure 19. Scatter plot of q against do for microlensing events with de- 
tectable binary signatures. Caustic crossing events are plotted with red 
squares, and smooth events with blue circles. Events classified as orbital 
motion events are plotted with larger, darker filled points, and those clas- 
sified as static with smaller, lighter, open points. The black lines show the 
positions of the caustic topology boundaries. 
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Figure 20. As Figure[l9]but showing |u u | plotted against do- 



events, and is not in itself evidence of two classes of orbital mo- 
tion events, but is instead a result of different caustic sizes in the 
different caustic topologies. 

Another feature of the plot is the clustering of caustic cross- 
ing orbital motion events near the boundary of the close and in- 
termediate topologies. It is close to the topology boundaries that 
the changes in projected separation cause the largest changes in the 
caustics. It is however difficult to attribute this clustering to faster 
caustic motions due to separational changes, as orbital velocity is 
inversely correlated with do, and so there should be more orbital 
motion events at smaller values of do in any case. In support of the 
existence of a separational class, there is a hint of clustering against 
the resonant-wide boundary. However, the caustic size peaks at 
both topology boundaries, as the single resonant caustic stretches 
before splitting apart into central and secondary caustics, possibly 
meaning that simply the increased size of the caustics causes the 
increased density of detections. 

Figure |20] plots the impact parameter against do and is very 
useful in separating different kinds of binary event, especially for 
planetary lenses. The events follow a distinctive pattern, with a 
large clump of events centred at |un| ~ and log do ~ which 
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consists of high-magnification events that encounter the central or 
resonant caustic. At very small |wo| this clump extends over a sig- 
nificant range in do, but narrows as [«o| increases, to its narrowest 
point at |uo| ~ 0.3, corresponding to the maximum size of the re- 
gion affected by resonant caustics (or at larger |uo| for stellar bina- 
ries). As | mo | increases, the plot shows a distinctive 'V shape, with 
no binary signatures being detected for events with do ~ 0. This 
'V shape arises as, in events with larger |uq|, the source passes 
through regions of the magnification pattern that can only con- 
tain secondary caustics, and does not enter the regions containing 
central or resonant caustics, i.e. the binary features in lenses with 
do ~ 1 only occur in regions of the magnification pattern that the 
sources with large |uo| do not probe. 

The events which occur on the branch with large | uo J and large 
do are caused by wide topology lenses, and therefore involve only 
a single secondary caustic encounter. The rotation of these lenses 
is typically very slow, and over the short duration of the binary 
features (typically of order a day), the lens completes only a very 
small fraction of its orbit. This points towards separational changes 
being the dominant effect in the detection of orbital motion features 
in events on this branch, even with the enhancement of rotational 
velocity due to the solid body 'lever arm'. 

The events that occur on the branch with large |uo | and small 
do are largely smooth events, with the occasional caustic crossing 
event. The smooth events are likely caused by the source crossing 
the large cusp extensions that occur in close binary lenses, suggest- 
ing that they will belong to the rotational class of events. 

Unfortunately it is difficult to attribute the cause of any one 
grouping of orbital motion events in Figures [19] and [20] to either 
the rotational or the separational class, partly because both types 
of motion will affect each event to some extent. Despite this, it is 
possible to classify many individual events as either a separational 
or rotational event. Figures [2l] and [22] show example lightcurves 
of both classes of orbital motion event, rotational and separational, 
respectively. The plots show the lightcurve in the upper left panels, 
with simulated data in red, the true model in blue, the best fitting 
static binary model in green and the best fitting single lens model in 
black. Also shown are the residuals from the single lens model and 
the static binary model in the middle and lower left panels respec- 
tively. Shown in the right panel is a plot of the source trajectory, 
shown in black, and snapshots of the caustics at various times dur- 
ing the event, shown in different colours. The coloured points on 
the time axis of the lightcurve show the time at which the caustic 
snapshots occurred, and the coloured points on the source trajectory 
show the position of the source at these times. The source trajectory 
and caustics are shown in the frame of reference that rotates with 
the binary axis, with its origin at the centre of mass. In this frame, 
rotation of the lens causes the source trajectory to appear curved, 
and changes in lens separation cause the caustics to change shape 
and move. Note that in event f in Figure |2T] and events e and f in 
Figure [22] the lens orbits are eccentric, so that the source does not 
travel along the shown trajectory at a constant rate. 

Figure |2T] shows examples of separational events. In each ex- 
ample the source trajectory appears relatively straight, indicating 
that the lens rotates little; however, in each case the caustics move 
significantly. Events a, b, c and e all involve resonant caustic cross- 
ings, and conform well to the picture described by Gaudi (2009). 
Event d could be described as the encounter of two disjoint caustics, 
similar to the original description of the rotational class of events 
by Gaudi ( 2009 1, but other than the close topology, the event is re- 
markably similar to event e; the source trajectory is slightly curved, 
but it is clear that separational effects are dominant. At first glance, 



event f would clearly fit into the picture of disjoint caustic encoun- 
ters, but the source trajectory reveals that rotation plays only a mi- 
nor role. In this event, a static fit to just the features about t = to 
would suggest a close encounter with a large secondary caustic at 
t « 1.5*e, but instead changes in the binary's separation cause 
the source to not just encounter, but cross a now much smaller sec- 
ondary caustic at t ~ 2£e- 

In contrast to Figure [21] the source trajectories in Figure [22] 
show significant curvature. Event a fits the description of rotational 
events by |Gaudi| p009l exactly. The source first encounters a sec- 
ondary caustic, but the rotation of the lens causes the source to 
pass the opposite side of the central caustic. Rotation also pre- 
vents the source from crossing the magnification excess between 
central caustic and the other secondary caustic. During the en- 
tire event, separational changes cause only slight changes in the 
caustics. In event c the rotation is more extreme, but the caus- 
tics smaller. The binary features are therefore more subtle, being 
caused by small magnification excesses between the caustics; the 
secondary caustics being located at ~ (—3, ±4), and the central 
caustic at ~ (0, 0). The rotation of the lens causes the source to 
cross each excess more than once, and there are several minor de- 
viations visible in the residual between the static and true model 
of the event. Event d, while being caused by a wide lens, expected 
to rotate slowly, is clearly caused by rotation. During the event, 
there are virtually no separational changes, but the precision with 
which the secondary caustic crossing and cusp approach features 
constrain the source trajectory mean that the very slight rotation 
which brings the source closer to the central caustic is detectable. 
Events b and e both show strong signs of rotation in their source 
trajectories, but separational changes are also important. While we 
assign them to the rotational class of event, in reality they may bet- 
ter fit into a third, hybrid class. Event f also shows signs of both 
rotational and separational orbital motion effects, but we assign it 
to the rotational class because without rotation the second caustic 
crossing would be significantly shorter. 



6 DISCUSSION 

6.1 Limitations of the study 

The questions that we wanted to answer in this work were: what 
fraction of microlensing events observed by the next generation 
surveys will be affected by orbital motion and what type of events 
are the effects likely to be seen in? While we do not claim to have 
fully answered these questions, we do feel that this work represents 
an important step in that direction. The simulation of the photome- 
try is slightly optimistic, and does not include the effects of weather 
and the systematic differences in the site conditions and observing 
systems, distributed across the Globe, that would make up the net- 
work of telescopes needed for a continuous monitoring microlens- 
ing survey. The observing setup we simulated is in some respects 
more like a space based microlensing telescope than a ground based 
network. However, the photometric accuracy that we simulated is 
not too optimistic, and the differences between the static and orbit- 
ing simulations show that orbital motion plays a significant role in 
a significant fraction of microlensing events. 

As discussed in Section [3] our choice of models will not fully 
answer the question of how many microlensing events with orbital 
motion effects will be seen, however, they do provide a good order 
of magnitude estimate. The binary detection efficiencies we find 
assume that all stars have a companion, and so must be adjusted 
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Figure 21. Example lightcurves of simulated events affected by separational type orbital motion effects. In each subfigure, the left panels show the lightcurve, 
its residual with respect to the best fitting Paczynski model and its residual with respect to the best fitting static binary model, from top to bottom respectively. 
Simulated data is shown in red, the Paczynski model in black, the static binary model in green and the true model in blue. The right panel shows the caustics 
at various times, and the source trajectory in the frame of reference rotating with the projected binary axis. The source trajectory is plotted in black, and the 
caustics are colour coded according to the time. Coloured points on the lightcurve panel show the time at which the caustic was in the state shown, and the 
coloured points on the source trajectory show the position of the source at this time. The parameters of the microlensing events can be found in Tables 4 and 5 
in the online supplementary material. 
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Table 4. Microlensing parameters of example lightcurves in the paper. 
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tFor events with eccentric orbits, two values of inclination are quoted, representing inclinations about two orthogonal axes on the sky. The effect of 



this second inclination is absorbed into the source trajectory for circular orbits, and to first order can be reduced to the range 0° ^ i ^ 90°. 



accordingly to account for this. For example, current estimates sug- 
gest that only ~ 33 percent of stellar systems are binaries (e.g. |Lada| 
2006), so assuming that a next generation microlensing survey de- 
tects ~ 2000 events per year we can expect to see ~ 30 stellar 
binary microlensing events showing orbital motion signatures per 
year. However, the true rate may be higher as the mass ratio distri- 
bution that we use for stellar binaries is not realistic; the real distri- 
bution is likely to be peaked in the range 0.1 4 (-1 4 1 (e.g. Duquen- 
|noy & Mayor 1991 1. A similar calculation for planetary lenses, as- 
suming the fraction of stars hosting planets is ~ 0.5, yields a detec- 
tion rate of ~ 1.5 caustic crossing orbital motion events per year. 
Again, this estimate is affected significantly by our assumptions. 
Our mass ratio distribution is optimistic, with current microlensing 
results suggesting an inverse relation between planet frequency and 
mass ratio in the regions microlensing is sensitive to (Su mi et al.| 
|2010| |Gould et al.|2010) . This implies our estimate is optimistic, 
but we have also assumed there is only one planet per system. Many 
multiplanet systems have been discovered to date (e.g. Gaudi et al. 
[20081 |Fischer et al.|2008> , and they are thought to be common. The 
microlensing planet detection efficiency in multiplanet systems is 



increased, as the planets are spread over a range of semimajor axes. 
This will somewhat compensate for the overestimate due to the in- 
correct mass ratio distribution. 

The major limitation of this work is that finite source effects 
are not considered. The finite size of the source acts to smooth out 
the extreme magnification peaks as a source crosses a caustic, lim- 
iting the precision with which magnifications can be measured, and 
caustic crossings timed, and thus plays an important role in orbital 
motion detection. However, in most cases, the caustic entry times 
can still be timed accurately if the caustic crossing is monitored 
with high enough cadence. In some cases the effect may increase 
the detectability of orbital motion as the source will probe more 
of the magnification pattern, especially when a source travels ap- 
proximately parallel to and very close to the inside of a fold caus- 
tic, producing additional peaks between the caustic crossings. We 
cannot quantitatively estimate the effects that finite source size has 
on the orbital motion detection efficiency, but we do not believe it 
will significantly affect our order of magnitude estimates. Unfortu- 
nately including finite source sizes in the modelling of a microlens- 
ing event increases the required computation time by several orders 
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of magnitude, so the effect could not easily be included in the sim- 
ulations without significantly reducing the sample size. 



6.2 Comparison with observations 

While our simulations are more representative of future microlens- 
ing surveys, it is possible for us to compare the results of our sim- 
ulations with the results of the current microlensing observations. 
Current microlensing planet searches using the survey/follow-up 
strategy routinely achieve a cadence similar or better than that ex- 
pected for future high-cadence surveys for a small number of mi- 
crolensing events per year (e.g. Dong et al. 2009 1. We can therefore 
compare the detection efficiency of orbital motion in the events 
where planets are detected. At the time of writing, there were 10 
published detections of planets by microlensing (Bond et al. 2004 
Udalski et al.|2005||Bea~ulieu et al.|2006||Gould et al.|2006t|Gaudi] 
et al.|2008||Bennett et al.T2008l[Donget al.|2009||Sumi et al.|20T0[ 



^ . 1 1 ■ 1 1 ■ ■ 

Janczak et al. 2010 ), and of these, 7 had high cadence coverage of 

a significant proportion of the lightcurve. In two of these events the 

orbital mot ion of the planet was detected ( Gaudi et al. 2008] [Dong] 

et al.|2009 Irlleading us to estimate an orbital motion detection ef- 



ficiency of ~ 0.29+q' ]^ percent. This efficiency is larger than we 
find in our simulations. However, the orbital motion effects in the 
OGLE-2005-B LG-71 event are very subtle, and improv e the fit by 
Ax! < 200 l |Udalski et al.||2005| |Dong et al.|| 2009 if] meaning 
that it would not be classed as a detection in our simulations; this 
reduces the comparable detection efficiency estimate to CI^q'^, 
with which our estimates of the detection efficiency for planetary 
caustic crossing events are consistent. It should be noted that this 
figure could be biased as events showing orbital motion signatures 
will take significantly longer to analyze. Unfortunately a similar es- 
timate for binary star lenses is not so simple as they are usually not 
followed-up to the same degree that planetary events are, either in 
terms of observations or modelling. 

We have identified two different classes of orbital motion 
event so it is natural to try to classify the orbital motion events 
that have already been seen. The orbital motion detected in OGLE- 
2006-BLG-109 ( |Gaudi et al.||2008l [Bermett et al.|[20T0t was de- 
tected due to deformation of a resonant caustic, and the event can 
easily be assigned to the class with separational changes. OGLE- 
2005-BLG-71 ( |Udalski et al.|2005[|Dong et al.|2009| > is harder to 
classify, as the orbital motion effects observed were very subtle. 
The event suffers from the well known close-wide degeneracy ( Gri- 
est & Safizadeh 1998 , Dominik 1999 ), and rather strangely, for the 
close (d < 1) solution, separational changes are more prominent 
than rotational, and vice versa for the wide (d > 1) solution, where 
we might normally expect the opposite. We therefore do not assign 
the event to either class. Of the stellar binary lenses, MACHO-97- 
BLG-41 (Albrow et al. 2000) was mainly influenced by rotation, 
and was detected by two disjoint caustic crossings, so is classed as a 
rotational event. EROS-2000-BLG-5 ( |An et al. 12002} undoubtedly 
belongs to the separational class; the caustic structure was resonant 



5 While the orbital motion of the Jupiter analogue was not detected in the 
OGLE-2006-BLG-109 system, the planet itself would still have been de- 
tected in the absence of the Saturn analogue, so it contributes to the denom- 
inator of the detection efficiency, but not to the numerator. 

6 The overall reduction in x 2 between the two analyses was much less 
than 200 when the size of the data sets and differing degrees of freedom are 
accounted for. The full analysis by Dong et al. 1 2009 1 included higher order 
effects not included in the original Udalski et al. (2005| analysis, some of 
which had a much larger effect than orbital motion. 



with d close to d w , and changes in separation were measured with 
high significance, while rotational changes were consistent with 
zero. The final events, OGLE-2003-BLG-267 and OGLE-2003- 
BLG-291 (Jaroszynski et al. 2005 ) are not very well constrained, 
so we do not attempt to classify them. 

We finally suggest that the event OGLE-2002-BLG- 
069 ( [Kubas et al.|2005| is a strong candidate for showing rotational 
type orbital motion effects. The event was modelled successfully 
by |Kubas et al.| ( |2005| > without including orbital motion, with a 
close binary solution favoured physically and by the modelling. 
The event had a time-scale £e ~ 105 d and binary parameters 
d — 0.46 and q — 0.58. The lightcurve was very similar to event 
b shown in Figure [22] having a long, well covered central caustic 
crossing, with measurements of both caustic entry and exit. The 
physical lens parameters obtained from the modelling suggest lens 
masses of Mi = 0.51Mq and M2 = 0.30Mq, and a projected 
separation of ~ 1.7 AU, with a corresponding minimum period of 
T > 900 d. The baseline is relatively bright, at lb ~ 16.2, and so 
subtle magnification deviations could probably be constrained by 
the data, if they have been covered. 



6.3 Future prospects 

Interestingly, our results show that the orbital motion detection ef- 
ficiency depends only weakly on the mass ratio. In the case of 
planetary events, caustic crossing orbital motion detections occur 
preferentially in high to moderate magnification events (A > 5), 
while smooth orbital motion detections occur in all but high- 
magnification events. Our results therefore suggest that the strategy 
of targeting high-magnification events (Griest & Safizade h|1998| 
|Han & Kim||2003"| > should allow caustic crossing orbital motion 
events to be detected efficiently. However, the strong dependence of 
orbital motion detection efficiency on the event time-scale suggests 
that long time-scale events should also be routinely followed-up. 
While follow-up of these events requires a significant investment of 
resources from the follow-up teams, like high-magnification events 
they are relatively rare. For a given cadence, these events allow a 
better signal to noise detection of planetary deviations, and also al- 
low more time for the prediction of future features. Long time-scale 
events are also more likely to show parallax features, allowing con- 
straints to be placed on the lens mass. 

High cadence, continuous monitoring microlensing surveys 
will begin operating in the next few years. Already, the MOA-II 
survey (Hearnshaw et al. 2005, Sako et al. 2008) has been sur- 
veying a fraction of its total survey area with a cadence of ~ 10 
minutes for some time, and the OGLE-IV survey I Udalsk i|2009| > 
has begun operations this year, and should provide significant in- 
creases in cadence over OGLE-III. KMTnet, a uniform network of 
telescopes with near continuous coverage, and operating at a ca- 
dence of ~ 10 minutes should begin operating around 2013; this 
promises an almost order of magnitude increase in the detection 
rate of microlensing events, and a similar if not bigger increase 
in the detection rate of planets by microlensing. The uniform na- 
ture of the survey network will also make statistical analysis of the 
planets detected easier, greatly enhancing the work already done in 
this direction jSumi et al.|2010[ |Gould et al.|2010| >. The work we 
have presented shows that a significant fraction of the events will 
show signs of orbital motion, which will significantly complicate 
the interpretation of future planet detections. However, these com- 
plications can be used to provide valuable additional constraints on 
the lens. 

Often overlooked are binary star microlensing events. The 
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next generation surveys will detect many more binary star events 
than planetary events. A large number of these lenses will be lo- 
cated in the Galactic bulge and be composed of low mass stars, 
providing an opportunity to study the properties of the bulge binary 
star population. Our results show that a significant fraction of these 
events will show orbital motion signatures, and it is likely that in 
a significant number of these events it will be possible to measure 
the masses of the system. It should therefore be possible to mea- 
sure the statistics of a population that is difficult to reach by current 
spectroscopic and astrometric methods due to their low brightness 
and long periods. 



7 SUMMARY 

We have simulated the lightcurves of ~ 100, 000 microlensing 
events caused by stars orbited by a companion star or planet. By 
fitting simulated data with single lens and static binary models we 
have determined the fraction of these events where the binarity of 
the lens is detected, and we also estimate the fraction of these events 
where orbital motion is detected. For an observational set up that 
resembles a near future microlensing survey conducted by a global 
network of telescopes without intensive follow-up observations, we 
found that orbital motion was detected in ~ 5-10 percent of sim- 
ulated binary star microlensing events depending on the character- 
istics of the event. Similarly, the rate of detection of orbital motion 
in simulated microlensing events where a planet is detected was 
~ 1-5 percent. 

We investigated the effects of various event parameters on the 
fraction of events showing orbital motion. Orbital motion detection 
efficiency as a fraction of binary detections was found to depend 
only weakly on the mass ratio of the binary, but strongly on the 
event time-scale. We found that a significant number of microlens- 
ing events showing orbital motion can be classified into one of two 
classes: those where the dominant cause of orbital motion effects is 
either the separational motion of the binary due to either inclination 
or eccentricity, or the rotational motion of the binary. 



ACKNOWLEDGMENTS 

MP acknowledges the support of an STFC studentship and would 
like to thank Scott Gaudi for several interesting and helpful discus- 
sions, Martin Smith for providing the blending data and invaluable 
suggestions regarding the false positives correction, and the anony- 
mous referee for many suggestions that improved the paper. 



REFERENCES 

Albrow M. D., Beaulieu J. -P., Caldwell J. A. R., Dominik M., 
Gaudi B. S., Gould A., Greenhill J., Hill K., Kane S., Martin R., 
Menzies J., Naber R. M., et al., 2000, ApJ, 534, 894 

Alcock C, Allen W. H., Allsman R. A., Alves D., Axelrod T. S., 
Banks T. S., Beaulieu S. F, Becker A. C, Becker R. H., Bennett 
D. P., et al., 1997, ApJ, 491, 436 

Alcock C, Allsman R. A., Alves D., Axelrod T. S., Baines D., 
Becker A. C, Bennett D. P., Bourke A., Brakel A., Cook K. H., 
et al., 2000, ApJ, 541,270 

Alcock C, Allsman R. A., Alves D., Axelrod T. S., Bennett D. P., 
Cook K. H., Freeman K. C, Griest K., Guern J., Lehner M. J., 
Marshall S. L., Peterson B. A., Pratt M. R., Quinn P. J., Rodgers 
A. W, Stubbs C. W, Sutherland W, 1995, ApJ, 454, L125+ 



An J. H., Albrow M. D., Beaulieu J.-P, Caldwell J. A. R., DePoy 

D. L., Dominik M., Gaudi B. S., Gould A., Greenhill J., Hill 
K., Kane S., Martin R., Menzies J., Pogge R. W., Pollard K. R., 
Sackett P. D., Sahu K. C, Vermaak P., Watson R., Williams A., 
2002, ApJ, 572, 521 

An J. H, Gould A., 2001, ApJ, 563, Lll 1 

Beaulieu J.-P, Bennett D. P., Fouque P., Williams A., Dominik M., 
J0rgensen U. G., Kubas D., Cassan A., Coutures C, Greenhill J., 
et al., 2006, Nature, 439, 437 

Bennett D. P., Anderson J., Bond I. A., Udalski A., Gould A., 
2006, ApJ, 647, L171 

Bennett D. P., Bond I. A., Udalski A., Sumi T., Abe F, Fukui A., 
Furusawa K., Hearnshaw J. B., Holderness S., Itow Y., et al., 
2008, ApJ, 684, 663 

Bennett D. P., Rhie S. H, Nikolaev S., Gaudi B. S., Udalski A., 
Gould A., Christie G. W, Maoz D., Dong S., McCormick J., 
etal.,2010, ApJ, 713, 837 

Binney J., Tremaine S., 2008, Galactic Dynamics: Second Edi- 
tion. Princeton University Press 

Bond I. A., Udalski A., Jaroszyhski M., Rattenbury N. J., 
Paczyhski B., Soszynski I., Wyrzykowski L., Szymahski M. K., 
Kubiak M., Szewczyk O., et al, 2004, ApJ, 606, L155 

Cassan A., 2008, A&A, 491, 587 

Dominik M., 1998, A&A, 329, 361 

Dominik M., 1999, A&A, 349, 108 

Dong S., Bond I. A., Gould A., Kozlowski S., Miyake N., Gaudi 
B. S., Bennett D. P., Abe F, Gilmore A. C, Fukui A., et al., 2009, 
ApJ, 698, 1826 

Dong S., Gould A., Udalski A., Anderson J., Christie G. W., Gaudi 
B. S., The OGLE Collaboration Jaroszyfiski M., Kubiak M., 
Szymahski M. K, Pietrzyhski G., Soszynski I., et al., 2009, ApJ, 
695, 970 

Duquennoy A., Mayor M., 1991, A&A, 248, 485 

Erdl H, Schneider P., 1993, A&A, 268, 453 

Fischer D. A., Marcy G. W., Butler R. P., Vogt S. S., Laughlin G, 
Henry G. W., Abouav D., Peek K. M. G, Wright J. T, Johnson 
J. A., McCarthy C, Isaacson H, 2008, ApJ, 675, 790 

Gaudi B. S., , 2009, Oral presentation at 13 th Microlensing Work- 
shop, Paris 

Gaudi B. S., Bennett D. P., Udalski A., Gould A., Christie G. W, 
Maoz D., Dong S., McCormick J., Szymahski M. K., Tristram 
P. J., Nikolaev S., Paczyhski B., et al, 2008, Science, 319, 927 

Gould A., 1992, ApJ, 392, 442 

Gould A., Dong S., Gaudi B. S., Udalski A., Bond I. A., Greenhill 
J., Street R. A., Dominik M., Sumi T, Szymahski M. K, et al, 
2010, ApJ, 720, 1073 

Gould A., Loeb A., 1992, ApJ, 396, 104 

Gould A., Udalski A., An D., Bennett D. P., Zhou A.-Y, Dong S., 
Rattenbury N. J., Gaudi B. S., Yock P. C. M., Bond I. A., et al., 
2006, ApJ, 644, L37 

Griest K, Hu W., 1992, ApJ, 397, 362 

Griest K, Safizadeh N., 1998, ApJ, 500, 37 

Han C, 2006, ApJ, 638, 1080 

Han C, Gould A., 1995a, ApJ, 449, 521 

Han C, Gould A., 1995b, ApJ, 447, 53 

Han C, Kim Y, 2001, ApJ, 546, 975 

Hearnshaw J. B., Abe F, Bond I. A., Gilmore A. C, Itow Y, 
Kamiya K, Masuda K, Matsubara Y, Muraki Y, Okada C, 
et al., 2005, ArXiv Astrophysics e-prints 

Holtzman J. A., Watson A. M., Baum W. A., Grillmair C. J., Groth 

E. J., Light R. M., Lynds R., O'Neil Jr. E. J., 1998, AJ, 1 15, 1946 
Hwang K, Udalski A., Han C, Ryu Y, Bond I. A., Beaulieu J., 



Orbital Motion in Microlenses 21 



Dominik M., Home K., Gould A., Gaudi B. S., et al., 2010, ApJ, 
723, 797 

Ioka K., Nishi R., Kan-Ya Y., 1999, Progress of Theoretical 

Physics, 102, 983 
James R, Roos M., 1975, Computer Physics Communications, 10, 

343 

Janczak J., Fukui A., Dong S., Monard L. A. G., Kozlowski S., 
Gould A., Beaulieu J. P., Kubas D., Marquette J. B., Sumi T., 
et al., 2010, ApJ, 711,731 

Jaroszynski M., 2002, Acta Astronomica, 52, 39 

Jaroszynski M., Skowron J., Udalski A., Kubiak M., Szymanski 
M. K., Pietrzynski G., Soszynski I., Zebrun K., Szewczyk O., 
Wyrzykowski L., 2006, Acta Astronomica, 56, 307 

Jaroszynski M., Udalski A., Kubiak M., Szymanski M., Pietrzyn- 
ski G., Soszynski I., Zebrun K., Szewczyk O., Wyrzykowski L., 
2004, Acta Astronomica, 54, 103 

Jaroszynski M., Udalski A., Kubiak M., Szymanski M. K., 
Pietrzynski G., Soszynski I., Zebrun K, Szewczyk O., 
Wyrzykowski L., 2005, Acta Astronomica, 55, 159 

Kim S., Park B„ Lee C, Yuk I., Han C, O'Brien T„ Gould A., 
Lee J. W., Kim D., 2010, in Society of Photo-Optical Instrumen- 
tation Engineers (SPIE) Conference Series Vol. 7733 of Society 
of Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, Technical specifications of the KMTNet observation sys- 
tem 

Kubas D., Cassan A., Beaulieu J. P., Coutures C, Dominik M., 
Albrow M. D., Brillant S., Caldwell J. A. R., Dominis D., Dona- 
towicz J., et al., 2005, A&A, 435, 941 

Lada C. J., 2006, ApJ, 640, L63 

Mao S., Paczynski B., 1991, ApJ, 374, L37 

Mao S., Paczynski B., 1996, ApJ, 473, 57 

Nemiroff R. J., Wickramasinghe W. A. D. T, 1994, ApJ, 424, L21 

Newcombe R., 1998a, Statistics in Medicine, 17, 873 

Newcombe R., 1998b, Statistics in Medicine, 17, 857 

Paczynski B., 1986, ApJ, 304, 1 

Paczynski B., 1997, ArXiv Astrophysics e-prints 

Park B.-G., DePoy D. L., Gaudi B. S„ Gould A., Han C, Pogge 

R. W., Abe E, Bennett D. P., Bond I. A., Eguchi S., et al., 2004, 

ApJ, 609, 166 
Rahvar S., Dominik M., 2009, MNRAS, 392, 1 193 
Rattenbury N. J., Bond I. A., Skuljan J., Yock P. C. M., 2002, 

MNRAS, 335, 159 
Refsdal S„ 1966, MNRAS, 134, 315 

Ryu Y, Han C, Hwang K., Street R., Udalski A., Sumi T, Fukui 
A., Beaulieu J., Gould A., Dominik M., et al., 2010, ApJ, 723, 
81 

Sako T, Sekiguchi T, Sasaki M., Okajima K, Abe E, Bond I. A., 
Hearnshaw J. B., Itow Y, Kamiya K, Kilmartin P. M., Masuda 
K, Matsubara Y, Muraki Y, Rattenbury N. J., Sullivan D. J., 
Sumi T, Tristram P., Yanagisawa T, Yock P. C. M., 2008, Ex- 
perimental Astronomy, 22, 5 1 
Schechter P. L., Mateo M., Saha A., 1993, PASP, 105, 1342 
Skowron J., Jaroszynski M., Udalski A., Kubiak M., Szymanski 
M. K, Pietrzynski G., Soszynski I., Szewczyk O., Wyrzykowski 
L., Ulaczyk K, 2007, Acta Astronomica, 57, 281 
Smith M. C, Mao S., Paczynski B., 2003, MNRAS, 339, 925 
Smith M. C, Wozniak P., Mao S., Sumi T, 2007, MNRAS, 380, 
805 

Sumi T, 2004, MNRAS, 349, 193 

Sumi T, Bennett D. P., Bond I. A., Udalski A., Batista V., Do- 
minik M., Fouque P., Kubas D., Gould A., Macintosh B., et al., 
2010, ApJ, 710, 1641 



Sumi T, Wozniak P. R., Udalski A., Szymanski M., Kubiak 
M., Pietrzynski G., Soszynski I., Zebrun K, Szewczyk O., 
Wyrzykowski L., Paczynski B., 2006, ApJ, 636, 240 

Udalski A., 2003, Acta Astronomica, 53, 291 

Udalski A., 2009, in K. Z. Stanek ed., Astronomical Society of 
the Pacific Conference Series Vol. 403 of Astronomical Society 
of the Pacific Conference Series, The Optical Gravitational Lens- 
ing Experiment (OGLE): Bohdan's and Our Great Adventure, pp 
110-+ 

Udalski A., Jaroszynski M., Paczynski B., Kubiak M., Szymanski 
M. K, Soszynski I., Pietrzynski G., Ulaczyk K, Szewczyk O., 
Wyrzykowski L., et al„ 2005, ApJ, 628, L109 

Wilson E. B., 1927, Journal of the American Statistical Associa- 
tion, 22, 209 

Witt H. J., 1990, A&A, 236, 311 

Witt H. J., Mao S., 1994, ApJ, 430, 505 



